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BACKGROUND OF THE INVENTION 

This disclosure relates, in general, to a computer method. More particularly, this 
disclosure relates to a computer method for the analysis of communications networks that 
10 approximates the total cost of network bandwidth services on a Virtual Private Network. 
„ BACKGROUND OF THE INVENTION 

Cj The operating cost for a digital telephony network is relatively easy to determine when 

M the network operator implements the network at the physical layer. The total cost is the sum of 
yQ the cost of the hardware and the cost of the link. The cost of the hardware (e.g., switches, 
5 15 routers, or multiplexers) is a combination of the depreciated value of the hardware purchase price 
Li and the hardware maintenance charges. The cost of the link depends on whether the network 
==2 operator leases or owns the trunk lines. If the network operator leases the lines, the cost of the 
link is the sum of the monthly lease rate for the line. If the network operator owns the lines, the 
cost of the link is the depreciated purchase price of the line. 
20 A Virtual Private Network (VPN) shares resources, such as the facilities of a public 

network, with other virtual private networks. The VPN maintains privacy by encapsulating the 
VPN messages in packets or frame relay frames to separate the messages from those of other 
VPNs that share the same public network. Tunneling involves the process of using 
encapsulation to carry one protocol over another. Encrypting the messages in the packets can 
25 further enhance the privacy of the VPN messages. Additional discussion of virtual private 
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networks is provided in the book by David E. McDysan entitled "VPN Applications Guide", 
published by John Wiley & Sons, 2000. 

Replacement of the physical links with Virtual Circuits is a subsequent advance in 
technology that allows a telephony service provider to easily construct a network. A Permanent 
Virtual Circuit (PVC) uses a fixed logical channel to maintain a permanent association between 
data terminals. Since PVCs do not require a setup operation before sending data nor a 
disconnect operation after sending data, a new connection between the same users may route 
messages along a different path. The Virtual Circuits derive bandwidth from multiplexers using 
Time Division Multiplexing or by connections through a Frame Relay or Asynchronous Transfer 
Mode (ATM) switch. The latter case usually involves a published tariff, such as, the cost to 
maintain the physical network connection, set-up fees, teardown fees, or usage fees. There are 
various ways to optimize a Permanent Virtual Circuit (PVC) such as by sharing a single access 
line. The total cost to operate the network, however, is still relatively easy to determine. 

Recent technology advances allow a telephony service provider to construct a Virtual 
Private Network (VPN) that does not include any PVCs. Instead, each node in the network has a 
Committed Access Rate (CAR) for ingress and egress traffic. The network does not require a 
destination for the traffic, but the traffic must remain within the VPN. As long as the traffic 
obeys the flow rate limitations specified by the CAR, the node will maintain a certain quality of 
service, that is, loss and delay. 

Multi-Protocol Label Switching (MPLS) is an emerging Internet Engineering Task Force 
(IETF) standard for constructing a data packet transit network. An MPLS network moves data 
packets from an ingress node to an egress node and relies upon conventional routers using legacy 
Local Area Network techniques to deliver the data packets to the end systems. In the MPLS 
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network protocol, each incoming packet includes a label that has an appearance similar to a 
packet header. As the packets traverse the label switch path, each node along that path makes a 
forwarding decision for the packet based solely on the contents of the label. As each node 
forwards a packet, the network strips the existing label from the message and adds a new label 
that provides forwarding instructions for the next node on the label switch path. The MPLS 
standard, dated August 10, 2000, is available on the IETF Web site at 

http://www.ietf.org/internet-drafts/drafl-ietf-mpls-arch-07.txt. Additional discussion of MPLS is 
provided in the book by Bruce S. Davie et al., entitled "MPLS: Technology and Applications", 
published by Morgan Kaufmann Publishers, 2000. 

MPLS enables the creation of "virtual circuits", or tunnels, across an Internet Protocol 
(IP) network. These virtual circuits make is possible to use MPLS to provision VPN services 
such as traffic engineering. Using MPLS for a VPN provides traffic isolation, much like an 
Asynchronous Transfer Mode (ATM) or frame relay service. A description of a method for 
applying using MPLS to provision a VPN on an IP network is available as a Request For 
Comment (RFC) document on the IETF Web site at http://www.ietf.org/rfc/rfc2547.txt. 

In a MPLS VPN, there is a committed access rate (CAR) for both the inbound and 
outbound traffic from each site. The destination of this traffic is unspecified, but the traffic must 
remain within the VPN. The VPN customer contracts for VPN service within the limits 
specified by the CAR. The provider of the VPN on a public network must price this service to 
meet the competitive conditions of the market place and to recover the provider's investment in 
the infrastructure of the network. The problem to solve is to develop a way to establish a price 
for the bandwidth that each VPN customer uses, knowing only the access rates, not the 
destination of the traffic. 
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SUMMARY OF THE INVENTION 

The invention disclosed herein is a computer method for approximating the total cost of 
the network bandwidth services on a Virtual Private Network that employs Multi-Protocol Label 
Switching technology. The algorithm computes a greatest lower bound and a least upper bound 
for the network bandwidth services. The invention uses these bounds to approximate the total 
cost of the network bandwidth services. 

The process for computing a greatest lower bound for the bandwidth used in a VPN that 
uses MPLS, includes the steps of: 

(a) Creating an auxiliary graph of Start and Finish nodes and one node 
for each node in the VPN that is saturated with inbound flow; 

(b) Creating a directed edge between the nodes representing a least 
cost of adding flow from a node having flow left to place to a node having 
capacity to receive the flow (looping for all edges thus created); 

(c) Calculating a shortest path algorithm to find the path from the Start 
node to the Finish node that gives the optimal way of serving the 
previously unserved unit of flow; and 

(d) Outputting L as a greatest lower bound for the bandwidth used. 

The process for computing in a separate calculation a least upper bound for the 
bandwidth used in a VPN that uses MPLS, includes the steps of: 

(e) Creating an auxiliary graph consisting of the nodes Start, Finish, 
and one node for each node in the VPN that is saturated with inbound 
flow; 
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(f) Creating a directed edge between the nodes representing a greatest 
cost of adding flow from a node having flow left to place to a node having 
capacity to receive the flow (looping for all edges thus created); 

(g) Calculating a longest path algorithm to find the path from Start to 
Finish, giving a worst way of serving the previously unserved unit of flow; 
and 

(h) Outputting U as a least upper bound for the bandwidth used. 

The process for combining the upper bound and the lower bound to produce an estimate 
for the cost of the bandwidth, includes the step of: 

(i) Setting a TotaLcost «axU + (l-a)xL 

where a is a measure of risk of losing money by carrying a given customer's VPN. 
DESCRIPTION OF THE FIGURES 

Figure 1 is a diagram of a complete weighted graph of a six node MPLS-based VPN 
network. 

Figure 2 is a flow diagram of the steps in determining a greatest lower bound for the 
bandwidth used in an MPLS-based VPN network. 

Figure 2(a) shows the start node "S" and the finish node "Fin" of the auxiliary graph for 
the first iteration in determining L as a greatest lower bound for the bandwidth used. 

Figure 2(b) shows the start node "S", the intermediate node "D" and the finish node "Fin" 
of the auxiliary graph for the second iteration in determining L as a greatest lower bound for the 
bandwidth used. 
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Figure 2(c) shows the start node "S", intermediate nodes "D" and "E", and the finish node 
"Fin" of the auxiliary graph for the third iteration in determining L as a greatest lower bound for 
the bandwidth used. 

Figure 3 depicts the steps in determining a least upper bound for the bandwidth used in an 
5 MPLS-based VPN network. 

Figure 3(a) shows the start node "S" and the finish node "Fin" of the auxiliary graph for 
the first iteration in determining U as a least upper bound for the bandwidth used. 

Figure 3(b) shows the start node "S" and the finish node "Fin" of the auxiliary graph for 
the second iteration in determining U as a least upper bound for the bandwidth used. 
44) Figure 3(c) shows the start node "S", intermediate node "A", and the finish node "Fin" of 

q the auxiliary graph for the third iteration in determining U as a least upper bound for the 
yl bandwidth used. 

ry Figure 3(d) shows the fourth iteration where flow is added from C to B. 

H= Figure 3(e) shows the fifth iteration in which an example is seen of moving flow. 

^5 Figure 4 shows the results as a plot of the number of networks vs. the cost per kilobit for 

'% fifty networks of five nodes each. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The invention is a computer method for approximating the total cost of the network 
bandwidth services on a Virtual Private Network (VPN) that employs Multi-Protocol Label 

20 Switching (MPLS) technology. A computer system for carrying out the method of the invention 
is a standard general-purpose data processor that includes a random access memory to store the 
program embodiment of the invention and a central processor to execute the instructions in the 
program embodiment. The computer system accesses databases that store the characteristics of 
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the MPLS-based VPN networks analyzed by the computer method invention. The computer 
system also connects to a billing network to report the total cost of the network bandwidth 
services on the MPLS-based VPN networks analyzed by the computer method invention. The 
billing network can connect to a banking system that debits payments from a VPN customer's 
network services account. 

Figure 1 is a diagram of a complete weighted graph that is an example of an MPLS-based 
VPN with six nodes that the computer method invention will analyze. Each of the network 
nodes labeled A, B, C, D, E, F, and G, has a possible communications path to every other node 
of the network. Each node is characterized by two values, an ingress flow rate capacity I, and an 
egress flow rate capacity O, which are expressed in units such as kilobits per second (Kbps). 
Node A, for example, has an ingress flow rate capacity I A = 20000 Kbps, and an egress flow rate 
capacity Oa = 20000 Kbps. Node B, for example, has an ingress flow rate capacity Ib = 20000 
Kbps, and an egress flow rate capacity Ob = 20000 Kbps. The path between node A and node B, 
which is called an edge, has a data transfer cost cab associated with it, which is expressed in units 
of dollars per Kbps. The data transfer cost c A b associated with the edge from node A to node B 
in this example is Cab = $3.50 per Kbps. Table 1 gives ingress flow rate capacity I and an egress 
flow rate capacity O for the network nodes labeled A, B, C, D, E, F, and G. Table 2 gives the 
data transfer cost for each of the possible communications paths between network nodes A, B, C, 
D, E, F, and G. 

Graph theory is well known in the art, and is described for example in the book by Frank 
Harary entitled "Graph Theory", published by Addison- Wesley, 1969. Using graph theory 
notation, a graph G is a pair 

G = (V,E) 9 
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where Vis a, finite set of nodes or vertices and E is a finite set that has as elements subsets of V 

of cardinality two called edges. The vertices of Fare usually called v h v 2 , By way of 

example, Figure 1 depicts a network that represents the graph 

f 

Weighted graphs are also known in the art. A graph G is a weighted graph if there is a 
real number associated with the edges in the graph. The weight of the edge from vertex vi to 
vertex v 2 is denoted as w ]2 . By way of example, Figures 2(a), 2(b), and 2(c) depict the network 
as a weighted graph. 

It is further known in the art that a graph G with n vertices in which each vertex is 
connected to each of the remaining vertices with one edge between each pair of vertices is a 
complete graph. 

The model of the network for the preferred embodiment is a complete weighted graph 

G = {V 9 (E, w E9 c E )), 

where Vis a finite set of network nodes, E is a finite set of edges (i.e., links) that connect pairs of 
nodes in set V 9 we represents the Administrative Cost associated with each link in set E, and ce is 
the cost per bit associated with each element of set E. If the network contains any parallel links, 
the model will collapse the parallel links into a single link with the associated costs for that 
single link equal to the sum of the cost of the parallel links. 

The Administrative Cost for a network link is used to make routing decisions and is 
usually set by the network operator based upon the recommendations of the hardware 
manufacturer. The cost per bit for a network link is a function of the operating point of the link 
and the recurring costs for the link. For example, the cost per bit for a Tl link (maximum 
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transfer capacity of 1,536,000 bits per second) between vertices vi and v 2 that is 65% utilized, 
only accommodates traffic flow in one direction, and costs $2200 per month is 

$2200 

° n ~ 1,536,0006^x0.65 ' 

If that same link can accommodate traffic in both directions, the cost per bit is 

$2200x2 
C]1 ~ 1,536,0006^x0.65 * 

The network for the preferred embodiment is a Virtual Private Network with a finite set of 
ingress rates / and egress rates O for each network node. At any particular time t 9 the traffic T 
being sent from node vpni to vprij is denoted as T.j(t). In dealing with stationary traffic, the 
traffic is denoted as 7] . . For any given pattern of the T Uj that satisfies the inequalities 

j 

The total cost of the flow is represented as: 

Total Cost = £X xC /j 

where, C. . is the sum of the cost per bit for each link on the shortest path from vpni to vprij as 

determined by the Administrative Cost for each edge. 

In the preferred embodiment, approximating the Total Cost requires calculating the least 
upper bound [/ m/w and greatest lower bound L max such that 

Total Cost « (a x C/ min )+((l - a)x L m J 
where the choice of a is a business decision for the telephony service provider. For larger 
values of a, the business decreases the risk of losing money by carrying a given customer's VPN 
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services and increases the risk that they will lose the contract to an aggressive competitor who 
undercuts their bid for the VPN services. For smaller values of a, the more the risks are 
reversed. The business must determine the value of a by instrumenting the VPNs and 
determining the actual value of a for a particular segment of the customer community. For 

example, if a telephony service provider determines that a = 0.4 with a variance of 0.1 for a 
given customer segment, a business decision to offer a new customer who is a member of that 
segment a price quote based upon a = 0.4 would be based on a reasonable confidence that any 
error in the calculation will be modest. 

In the network model of the preferred embodiment, the calculation of the least upper 
bound Umin and greatest lower bound L max is determined by solving, respectively, the longest and 
shortest path through the network. 

The computer method invention for approximating the total cost of the network 
bandwidth services on a VPN that employs MPLS, begins with determining a greatest lower 
bound for the bandwidth used, as shown in the flow diagram of Figure 2. 

The general process of computing a greatest lower bound for the bandwidth used in an 
MPLS-based VPN, as shown in the flow diagram of Figure 2, comprises the steps of: 

(A) Creating an auxiliary graph of Start and Finish nodes and one node 
for each node in the VPN that is saturated with inbound flow; 

(B) Creating a directed edge between the nodes representing a least 
cost of adding flow from a node which has flow left to place, to a node 
which has capacity to receive the flow (looping for all edges thus created); 

(C) Calculating a shortest path algorithm to find the path from Start 
and Finish nodes giving an optimal way of serving the nodes that have 
previously unserved units of flow; and 

(D) Outputting L as a greatest lower bound for the bandwidth used. 
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The details of calculating the Greatest Lower Bound (Lmax), are as follows. 

The calculation of the greatest lower bound is performed with a digital computer 
programmed to execute the following steps: 

1 . Initialize the Lower Bound L max to 0. 

2. Initialize the ingress flow and egress flow for each node in the VPN to 0. 

[a] Create an auxiliary graph that consists of a Start node, a Finish node, and each node in the 
VPN that has a saturated ingress flow. 

[b] Create a directed edge between every node in the auxiliary graph on the path from Start to 
Finish (looping for all edges thus created). 

[c] The directed edge from Start to Finish has the label (k, c\a, I) and weight c# that represents 
the least cost of adding flow from vptik (which has flow remaining to place) to vpni 
(which has the capacity to receive flow). 

[d] The directed edge from Start to saturated node vpnj has the label (k, cui, I) and weight cuj 
that represents the least cost of adding from vpn k (which has flow remaining to place) to 
the saturated node vpni (which has the capacity to receive flow). 

[e] The directed edge from saturated node vpnj to saturated node vpw/has the label (k, (cu- 
Cki)> J) an d weight (c^- cm) that represents the increased least cost of moving a unit of 
traffic from the saturated node vpni to the saturated node vpnj. 

[f] Run a Shortest Path algorithm to determine the path from Start to Finish that provides the 
optimal way of serving the previously unserved unit of flow. 

[g] Calculate L 1 as the sum of the cost for each edge in the path from Start to Finish. Add L ' to 
the current value of L max . • 

[h] Destroy the Auxiliary graph created in step [a] above. 
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[3] Repeat steps [2] until all of the demand in the system is allocated. 
[4] Output L max and terminate. 

The following example calculates the greatest lower bound in a VPN that has six nodes 
(A, B, C, D, E, and F). Table 1 describes the ingress and egress flow rates for each node. Table 
5 2 describes the cost matrix for transferring 1 kilo-bit of data traffic in one second between nodes 
in the VPN. 



Node 


Ingress Flow 


Egress Flow 




Rate 


Rate 


A 


20000 


20000 


B 


20000 


20000 


C 


20000 


20000 


D 


10000 


10000 


E 


10000 


10000 


F 


10000 


10000 



Table 1 : Ingress and Egress Flow Rates in the VPN. 
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Node 


A 


B 


C 


D 


E 


F 


A 


00 


^ so 


s so 


4 on 


6 00 

U. W 


R 00 


"D 
D 


1 SO 


CO 


0 00 
y.vjyj 


R 00 
o. w 


1 70 


8 00 


c 


s so 

J. JU 


y.\j\j 


CO 


7 00 


6 00 


7 SO 


D 


4.00 


8.00 


7.00 


oo 


2.00 


2.50 


E 


6.00 


3.70 


6.00 


2.00 


CO 


3.00 


F 


8.00 


8.00 


7.50 


2.50 


3.00 


oo 



Table 2: Data Transfer Costs ($/kbps) to Compute L 



After initializing the lower bound and the ingress flow and egress flow for each node in 
the VPN, the digital computer program creates the auxiliary graph depicted in Figure 2(a). 
Figure 2(a) shows the start node "S" and the finish node "Fin". The digital computer program 
creates a directed edge from node E (Start) to node D (Finish) and assigns a 10 kbps traffic flow 
rate to that edge. Since the shortest path is from node E to node D, 
$2 00 

V= — — x 10 kbps = $20.00 . The value of L after the first iteration through the algorithm is 
kbps 

$20.00. 

Since the first iteration saturated the ingress flow to node D } the digital computer 
program creates the auxiliary graph depicted in Figure 2(b) for the second iteration. Figure 2(b) 
shows the start node "S", the intermediate node "D" and the finish node "Fin". The digital 
computer program creates a directed edge from node D (Start) to node E (Finish) and assigns a 
10 kbps traffic flow rate to that edge. In addition, the digital computer program creates directed 
edges from node F to node D and from node E to node F. The path through node D has a cost of 
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$3.50 while the path, from node D to node E only has a cost of $2.00. Thus, 
$2 00 

L'= — — xlOkbps = $20.00 and the value of L after the second iteration through the algorithm 
kbps 

is $40.00. 

Since the first and second iterations saturated the ingress flow to node D and E, the 
digital computer program creates the auxiliary graph depicted in Figure 2(c) for the third 
iteration. Figure 2(c) shows the start node "S", intermediate nodes "D" and "E", and the finish 
node "Fin". This auxiliary graph consists of seven directed edges. Two of the edges, node D to 
node E and node E to node D, have weight oo because the flow into the one node originates from 
the other node and the disclosed network model does not allow for traffic flow from a node onto 
itself. The auxiliary graph depicted in Figure 2(c) has three shortest paths. At this point, the 
algorithm's course is implementation dependent. In the preferred embodiment, the algorithm 
chooses to add 20 Kbps of traffic flow from node B to node A. Thus, 
$3 50 

V- — : — x 20kbps = $70.00 and the value of L after the third iteration through the algorithm is 
kbps 

$110.00. 

The algorithm continues through a four additional iterations each time saturating the 
ingress flow into another node and adding another node to the auxiliary graph. In the fourth 
iteration, for example, the auxiliary graph consists of the nodes Start, D, E, A, and Finish. Table 
3 summarizes the results after each iteration. After completion of the iterations, the value of L 
is the optimal flow because each step in the iteration produced an optimal flow. 
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Iteration 


# of Nodes in 
the Auxiliary 
Graph 


Flow Rate 
(Directed Edge) 
Added 


L' 


L 


1 


2 


10K(DtoE) 


20 


20 


2 


3 


10K(EtoD) 


20 


40 


3 


4 


20K(BtoA) 


70 


110 


4 


5 


20K (A to B) 


70 


180 


5 


6 


-10K (E to D) 
10K (E to F) 
10K (FtoD) 


35 


215 


6 


7 


-10K (E to F) 
10K(EtoC) 
-lOK(DtoE) 
10K(DtoF) 
10K(CtoA) 


90 


305 


7 


7 


-10K(BtoA) 
lOK(BtoC) 
lOK(CtoA) 


110 


415 



Table 3: Summary of the Iterations in the Example for Computing the Lower Bound 
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The general method of the invention for approximating the total cost of the network 
bandwidth services on a VPN that employs MPLS, continues by determining a least upper bound 
for the bandwidth used. This is generally shown in Figure 3. 

The general process of computing a least upper bound for the bandwidth used in an 
MPLS-based VPN, as shown in Figure 3, comprises the steps of: 

(E) Creating an auxiliary graph consisting of the nodes Start, Finish, 
and one node for each node in the VPN that is saturated with inbound 
flow; 

(F) Creating a directed edge between the nodes representing a greatest 
cost of adding flow from a node which has flow left to place, to a node 
which has capacity to receive the flow (looping for all edges thus created); 

(G) Calculating a longest path algorithm to find the path from Start to Finish, giving a 
worst way of serving the nodes that have previously unserved units of flow; and 

(H) Outputting U as a least upper bound for the bandwidth used. 

Figure 3 also shows the remaining step of combining the upper bound and the lower 
bound to produce an estimate for the cost of the bandwidth, which is the following step of: 

(I) setting a Totalxost «axU + (l-a)xL 

where a is a measure of risk of losing money by carrying a given customer's VPN. 

The value of Totalxost is output by the computer to a billing network to report the total 
cost of the network bandwidth services on the MPLS-based VPN networks being analyzed by the 
computer method invention. The billing network can be connected to a banking system, for 
example, that debits payments from a VPN customer's network services account. 

The details of calculating the Least Upper Bound {U min ) are as follows: 

The calculation of the least upper bound is similar to the calculation of the greatest upper 

bound, however, the cost matrix is modified to set the cost of the diagonal terms to 0, not oo. 

Table 4 contains the cost matrix for calculating U. 
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Node 


A 


B 


C 


D 


E 


F 


A 


0 


3.50 


5.50 


4.00 


6.00 


8.00 


B 


3.50 


0 


9.00 


8.00 


3.70 


8.00 


C 


5.50 


9.00 


0 


7.00 


6.00 


7.50 


D 


4.00 


8.00 


7.00 


0 


2.00 


2.50 


E 


6.00 


3.70 


6.00 


2.00 


0 


3.00 


F 


8.00 


8.00 


7.50 


2.50 


3.00 


0 



Table 4: Data Transfer Costs ($/kbps) to Compute U 
O The calculation of U is performed with a digital computer programmed to execute the 

«j following steps: 

m 1. Initialize the Upper Bound U to 0. 

Hj5 2. Initialize the ingress flow and egress flow for each node in the VPN to 0. 

H [a] Create an auxiliary graph that consists of a Start node, a Finish node, and each node in the 

C : VPN that has a saturated ingress flow. 

U [b] Create a directed edge between every node in the auxiliary graph on the path from Start to 
Finish. 

10 [c] The directed edge from Start to Finish has the label (k, c k u I) and weight c# that represents 
the greatest cost of adding flow from vprik (which has flow remaining to place) to vpni 
(which has the capacity to receive flow), 
[d] The directed edge from Start to saturated node vpni has the label (k, cu, I) and weight c\a 
that represents the greatest cost of adding from vpnk (which has flow remaining to place) 

15 to the saturated node vpnj (which has the capacity to receive flow). 
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[e] The directed edge from saturated node vpnj to saturated node vnpj has the label (k, (cu - 
Cki), J) and weight (cm - cu) that represents the decreased greatest cost of moving a unit of 
traffic from the saturated node vpnj to the saturated node vpnj 

[f] Run a Longest Path algorithm to determine the path from Start to Finish that provides the 
5 least optimal way of serving the previously unserved unit of flow. 

[g] Calculate U' as the sum of the cost for each edge in the path from Start to Finish. Add U f to 
the current value of C/ m/ „. 

[h] Destroy the Auxiliary graph created in step [a] above. 

3. Repeat steps 2 until all of the demand in the system is allocated. 
Q0 4. Output U m i n and terminate. 

Ln The computer method invention for approximating the total cost of the network 

L j; bandwidth services on a VPN that employs MPLS, is completed by combining the upper bound 
1~ and the lower bound to produce an estimate for the cost of the bandwidth. This is done by the 
H5 step of: 

SJ setting a Totalxost «axU + (l-a)xL 

y where a is a measure of risk of losing money by carrying a given customer's VPN. 

20 

This algorithm is demonstrated by tracing through some of the steps. 

In the first auxiliary graph of Figure 3(a) the edge is labeled from S to Fin with a label of 
(F, 80, A). Notice that since the cost of going to B is the same, it is possible to also use the label 
(F, 80, B). At the end of this step, U= 80. The second auxiliary graph of Figure 3(b) also has the 
25 same node set since A can terminate 20 Kpbs of flow. The edge is given the label (F, 60, Q) and 
after this step, U= 140. 
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In the third step, since A is now saturated, the graph of Figure 3(c) has three nodes. The 
cost of adding the traffic from node D to node A is 40. The longest path is 80, thus the traffic is 
added from D to B, setting U to 220. In the fourth step shown in Figure 3(d), add 10 Kpbs of 
flow is added from C to B. In step 5 of Figure 3(e), the first example is seen of moving flow. 
5 Rather than adding 10 Kbps of flow from C to F at a cost of 75, it is better to first move the flow 
from C to B and then move the flow off of B onto C. As with the previous algorithm, not all of 
the stages are illustrated for the execution, but the remaining ones are left to the reader. 

Shortest Path (SP) algorithms such as Dijkstra are well known (see the Harary book, cited 
above), but Longest Path (LP) algorithms are not as well known. There is a good reason for this. 
CjO In general, LP problems are of exponential complexity. What saves this particular case is the 
^ special structure of the graph. It is possible to solve the LP problem for an auxiliary graph of 
1 2 node nbn — 2 executions of Dijkstras algorithm. Since each Dijkstra run is of complexity 
nj 0(n 2 ), the total complexity of finding the longest path from S to Fin isO(n 3 ). 

u Numeric Results 

□ 5 The previously introduced algorithms can be used to examine various simple network 

problems and give some insight into cost estimates. The model of Squareworld is used, which 
was introduced in the book by R. Cahn entitled "Wide Area Network Design: Concepts and 
Tools for Optimization", published by Morgan Kaufinann, 1998. Squareworld is a 1000 mile by 
1000 mile region in the United States. This lends itself to simple generation of network 

20 problems. 

Begin by assuming that the bandwidth cost is a simple linear function. I Kbps of 
bandwidth will cost 0.50 plus 0.10 per 100 km. Generate traffic at the nodes using a profile 
shown below: 
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2 /* 2 types of nodal traffic */ 

1 0000 1 0000 0.5 /* 50% of the nodes have 1 0 Kbps in and 1 0 Kbps out */ 

20000 20000 0.5 /* 50% of the nodes have 20 Kbps in and 20 Kbps out */ 

There are two types of traffic since these cases were used to debug the code and with one 
5 type of traffic the path through the code becomes simple matching of nearest and furthest 

neighbors. The results for 50 test runs are shown in Figure 4, which is a plot of the number of 
networks vs. the cost per kilobit for 50 networks of five nodes. 

The average cost runs anywhere from $0.835/Kbps to $1.665/Kbps showing that the ratio 
of the most expensive case is about twice the least cost case. 
EfO It is interesting to see what happens as the number of nodes increases. In an example 

^ where 50 experiments are run with 10 nodes in Squareworld., the lower bound is decreased from 
I n the previous case. Since it is known from the tariff structure that the cost for I Kbps is greater 
fu than or equal to $0.50, it is seen that the lower bounds for the individual flows are still well away 
from that bound. With 10 nodes, the average cost per Kbps varies between $1 .06 and $1 .48; thus 
45 with more sites the variance is decreased. 

u Extensions to MPLS for International VPNs 

Multi-protocol label switching (MPLS)-based VPNs do not adapt themselves particularly 
well to international use. Transoceanic telecommunications have cost between one and two 
orders of magnitude more than national telecommunications. To understand what this means for 
20 MPLS-based systems, two simple examples are reviewed. In the first model, it is assumed that 
there are 2n nodes that are in one country. It is assumed that cy = 1 and that // = 0i = 10 Kbps. In 
the second model it is assumed that nodes 1,. . ., n are in country A and that nodes n + 1,..., n are 
in country Z. It is assumed that eg = 1 if the nodes are in the same country and Cy = 20 if they are 
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in different countries. Again, /, = O l ? = 10 Kbps. Then if n is even, it can be can computed that 
L M0 dei 1 = U M0 dei 1 = 20/i. While L Mo dei2 = 20w, U M0 dei2 = 20n. Thus for any network with a large 
number of nodes on both sides of an ocean, the formula of 

Cost= a xL + (l-a)xU 

5 estimates the cost to be very large. 

The solution to this problem is to use the idea of a saturated cut. If the VPN was not 
provisioned with MPLS but rather provided as a set of PVCs and the sole connection between 
the United States and Europe was a 30 Kbps PVC, then the flow across the Atlantic would be 
limited by that link. The link is, in effect, a one-edge cut which disconnects the graph. If for 
3p reasons of reliability two-way connectivity is needed, replace the single PVC with a pair of 
□ PVCs, each with a committed information rate (CIR) of 15 Kbps. Then the two links still have a 

capacity of 30 Kbps and provide the same limit on the transatlantic flow. 
^ To offer cost effective multi-national MPLS-based networks, it will be necessary to 

[? modify MPLS to provide cut limits. The network will need, in effect, to not allow the flow 
3i5 between two sets of nodes to exceed a limit. Without this architectural addition, the cost (and 
O thus the price) of such MPLS-based VPNs will be prohibitive to the customer or risky to the 
provider. 

It is easy to modify the algorithms for computing L and U to take a cut constraint into 
account. This will result in a radical decrease in U and usually will involve no change to L at all. 
20 The algorithm to calculate U would then be: 

1. Sett/=0. 

2. Let Si and S2 be the two sets of nodes created by the cut. 

3. Loop until the cut is saturated. 
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(a) Find (s/, s 2 ) S\ x 52 such that C s i, S 2 is maximal, s/ can source 
additional flow and 82 has additional sink capacity. 

(b) U+ = c sIt s2 x FL 

(c) O sI -=FLI s2 - = FL 

4. End of loop. 

5. Compute £/;, the upper bound for the problem restricted to Si. 

6. Compute U 2 , the upper bound for the problem restricted to S2. 

7. U+ = U 1 + U 2 . 

The algorithm assigns the cut flow and then computes an upper bound for each side of the 
cut. While this does not produce a least upper bound, it is adequate for the purposes of estimating 
the tariff. 

Although a specific embodiment of the invention has been disclosed, it will be apparent 
to those having skill in the art that changes may be made to the specific embodiment without 
departing from the spirit and the scope of the invention. 
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